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Abstract 

> ■ 
m 

■ In this paper we study the statistical mechanics of the instanton hquid in 

O . QCD. After introducing the partition function as well as the gauge field and 

0^ 



quark induced interactions between instantons we describe a method to calcu- 
late the free energy of the instanton system. We use this method to determine 



Oh! 

^ • the equilibrium density and the equation of state from numerical simulations 



of the instanton ensemble in QCD for various numbers of flavors. We find that 

there is a critical number of flavors above which chiral symmetry is restored 

in the groundstate. In the physical case of two light and one intermediate 

mass flavor the system undergoes a chiral phase transition at T ~ 140 MeV. 

We show that the mechanism for this transition is a rearrangement of the 

instanton liquid, going from a disordered, random, phase at low temperatures 

to a strongly correlated, molecular, phase at high temperature. We also study 

the behavior of mesonic susceptibilities near the phase transition. 
PACS numbers: 12.38.Lg, ll.30.Rd, 05.70.-a 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Understanding the vacuum structure of gauge theories hke QCD is one of the main 
problems in quantum field theory today. It also provides the theoretical foundation for 
hadronic models and hadronic phenomenology from the underlying field theory of the strong 
interaction, Quantum Chromodynamics. There are a number of indications that instantons, 
classical tunneling trajectories in imaginary (euclidean) time are an important ingredient of 
the QCD vacuum. 

Soon after the discovery of instantons 20 years ago [T|], it became clear that instantons 
may provide at least a qualitative understanding of many features of the QCD vacuum. 
Instantons solve the U{1)a problem they give a mechanism for chiral symmetry breaking 
0, contribute to the gluon condensates and lead to a non-perturbative vacuum energy 0,^. 

The development of a quantitative theory based on these ideas took much longer. The 
instanton liquid model was originally suggested by Shuryak in 1982 0, guided mainly by 
phenomenological considerations. Later Diakonov and Petrov developed an analytical ap- 
proach based on the variational method (see also ||^). The first numerical simulations 
of the "instanton liquid" were reported in [p|,p!0|. During the past two years we have shown 



TT|-|T3| that the "random instanton liquid model" (RILM) provides a successful description 



of a large number of hadronic correlation functions, including mesons and baryons made of 
light quarks, heavy-light systems and glueballs. These correlators not only give reasonable 
values for the corresponding resonance masses and coupling constants, but they also com- 
pare well with point-to-point correlation functions extracted from phenomenology |]14| or 
measured on the lattice [|l^ . 

Following these developments, several recent lattice simulations have focused on the role 
of instantons in the QCD vacuum, both at zero and finite temperature. Using a method 
called "cooling", one can relax any given gauge field configuration to the closest classical 
component of the QCD vacuum. The resulting configurations were known to be of multi- 
instanton type |]T6|, but the more recent work by Chu et al. has provided quantitative 
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measurements of the parameters of the instanton hquid, as well as detailed studies of the 
dynamical effects of instantons. These authors conclude that the instanton density in the 
quenched theory (without dynamical fermions) at zero temperature is n ^ (1.3-1.6) fm~^ 
while the average size is about p ~ 0.35 fm. These numbers confirm the key parameters 
n = lfm~^ and p = 1/3 fm of the instanton liquid model mentioned above. In addition 
to that, Chu et al. studied correlation functions in the cooled configurations, finding that 
they hardly change from the original, fully quantum configurations. This would imply 
that instanton effects dominate over perturbative and confinement forces in determining 
the structure of low-lying hadronic states. A lattice measurement of the instanton size 



distribution for pure SU{2) gauge theory was performed in [jT8[, where also an attempt was 
made to study correlations between instantons. The size distribution was compared to the 
predictions of the interacting instanton liquid model in [0. 

In this paper we want to report a detailed study of the statistical mechanics of the 
interacting instanton liquid (IILM), both at zero and finite temperature. The purpose of 
this study is twofold. First, we want to give a much more consistent treatment of the model 
at zero temperature. Second, we explore the model at non-zero temperature. In this paper, 
we will limit ourselves to the analysis of bulk properties, like the energy density, the chiral 
condensate and mesonic susceptibilities. A detailed study of hadronic correlation functions 
will be presented in a forthcoming publication 

At zero temperature, we want to study the importance of instanton interactions and 
determine those features of the interacting ensemble that differ from the simplest, random 
ensemble (RILM). We also want to determine the relation between the instanton interaction 
and global parameters (like the density and average size) of the instanton liquid. Further- 
more, we want to quantify the role of instanton interactions in producing correlations among 
the instantons. This is of interest, since despite the success of the random model in the de- 



scription of most hadronic correlation functions [|TT|, it fails in channels where the single 
instanton interaction is strongly repulsive (like the t]' and 6 meson channels). As shown in 
2l| , [T3| , ^ , the correlations among instantons caused by their classical and fermion induced 



interactions lead to a correct description of topological charge screening and the 77' channel. 

The role of correlations among the instantons is particularly important with regard to 
the nature of the chiral phase transition. Originally, it was believed that chiral 



symmetry is restored due to the rapid disappearance of instantons at high temperature [25 



Ilgenfritz and Shuryak realized PB[ that instantons can be present even above the chiral 



phase transition, as strongly correlated instanton-antiinstanton molecules (see also [^]). In 
this case, the transition is determined by the phase equilibrium between the low temperature 
"liquid" and high temperature "molecular" phase. Shuryak and Velkovsky later argued that 
instanton suppression is essentially a plasma effect and should not be present below the 
phase transition p8[. This means that the transition is driven by the formation of molecules 
|l3T| , |32[| rather than by the suppression of individual instantons. There is some support for 
this scenario from lattice simulations performed by Chu and Schramm |2^. Extending the 



cooling method to finite temperature, they find that the instanton density is essentially 
independent of temperature below the phase transition, while it is exponentially suppressed 
above the transition temperature. First evidence for the presence of molecules near was 
reported in [BD . 



So far, the transition has only been studied using the schematic "cocktail" model in- 
troduced by Ilgenfritz and Shuryak |3l|. In this model the instanton liquid consists of two 
components, a random and a molecular one. The free energy is determined separately for 
the two components and their concentrations is then determined by minimizing the total 
free energy. The chiral phase transition occurs when the concentration of random instantons 
is zero. This approach predicts the presence of a substantial number of instantons even if 
T > Tc, causing new non-perturbative effects in the plasma phase. One such effect, studied 
p2| , is the "molecule-induced" effective interaction between quarks, leading to a spectrum 



m 



of spacelike screening masses consistent with lattice data. Another possible consequence. 



the survival of certain hadronic modes above the phase transition was studied in p3 |. 

In the present work we significantly improve on the schematic model used in these works, 
by doing a complete calculation in the interacting ensemble. In this way, many approxima- 



tions are relaxed and all possible correlations among instantons (not just polarized instanton- 
antiinstanton pairs) are included. The paper is organized as follows. In section 2 we intro- 
duce the partition function of the instanton liquid and specify the gauge field and fermion 
induced interactions between instantons. Detailed parameterizations of these interactions 
can be found in the appendices. In section 3 we describe the method which is used to calcu- 
late the partition function. In section 4 we use this method to study the instanton ensemble 
at zero temperature. In sections 5 and 6 we generalize the method to finite temperature 
and study the nature of the chiral phase transition in the instanton liquid. In section 7, we 
study some phenomena associated with the transition, in particular the Dirac spectrum and 
the mesonic susceptibilities. 

II. THE PARTITION FUNCTION OF THE INSTANTON LIQUID 

The euclidean partition function of QCD is given by 



where the gauge field action is given by S[Af^] = \J d'^xTi^F^^Ff^^) and the Dirac operator 

is defined by Dip = 7m(^m ~" i^n)"^- 

The main assumption underlying the instanton model is that the full partition function 
can be approximated by a partition sum in which the relevant gauge configurations are 
superpositions of instantons and antiinstantons. In this partition function the integration 
extends over the collective coordinates associated with instantons and instantons 



Here dQi = dUi d^Zi dpi is the measure in the space of collective coordinates, color orientation, 
position and size, associated with single instantons. For the gauge group SU (3), there are 12 
collective coordinates per instanton. Fluctuations around individual instantons are included 
in gaussian approximation. This gives the semi-classical instanton amplitude, originally 
calculated by 't Hooft 0. To two- loop accuracy, it reads 




/ 



(1) 




/ 
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d{p) = C^^p-'MPr""' exp (^-P,ip) + (2iV, -Yb^YbMp) 1°S(/51(P))) ' (3) 

4.6exp(-1.86iV,) 
7r2(Ar,- l)!(Ar,-2)! ' ^ 

where Pi{p) and P2{p) are the one and two loop beta functions 

A(p) = -Mog(pA), f3,{p) = f3,{p) + ^log(^/5i(p)), (5) 

with the one and two loop coefficients 

b = -N^- -Nf, b' = ^N'^ - ^-^N^Nf + ^. (6) 
3 3-^ 3 ^ 3 ^ A^c ^ ^ 

The coefficient C^^ was calculated in a Pauli-Vilars renormalization scheme, and A is the 
corresponding scale parameter. The classical action 5*0 = ^ is included in the semi- 
classical amplitude (|^). The classical interaction between instantons is denoted by Sint- We 
approximate this interaction by a pure two-body interaction S^t = | S/^j Sint{^ij) which 
only depends on the relative coordinates of the two instantons. The importance of genuine 
three body effects in the classical interaction between instantons was studied in , with the 
conclusion that this contribution is negligible as long as the density is not extremely large. 

The two-body interaction Sinti^ij) = S[Af^{Qjj)] — 2So is calculated classically, by 
inserting the two-instanton gauge potential A^{Qjj) into the action. There are no exact 
instanton-antiinstanton solutions to the classical Yang-Mills equations, so in practice one 
has to use an ansatz for the gauge potential. The resulting interaction will then depend on 
the details of the trial function. Various trial functions have been used in the literature: (i) 
the sum ansatz 0], (ii) the ratio ansatz (iii) the conformaly invariant Yung ansatz 



and (iv) the exact streamline solution |^. The latter is characterized by the fact that the 
action is minimized in all directions except along the collective coordinate describing the 
separation between the two instantons. In this sense, the streamline solution is the optimal 
classical instanton-antiinstanton configuration. 

In order to discuss the properties of the classical interaction between instantons, let us 
introduce the four vector = ^tr(f//t/AT+) where Uj^a are the orientation matrices of the 
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instanton and antiinstanton and r+ is the 2x2 matrix (r, —i). For the gauge group SU{2), 
is a real unit vector whereas for SU{3) it is a complex vector with \u\^ < 1. In any case, 
we can define an angle 6 by 

cos 6* = ' , , ' , (7) 

\u\ 

where R = zj — za is the vector connecting the centers of the two instantons. For all the trial 
functions mentioned above, the large distance part of the instanton-antiinstanton interaction 
is given by 

S.nt = P,{p)-^\u\'{l-4cos'e), (8) 

which is the dipole-dipole interaction originally discussed by Callan, Dashen and Gross 0. 
The interaction is given in units of the single instanton action f3i{p). The argument of the 
beta function is not exactly determined without a calculation of the fluctuations around the 
two-instanton configuration. In practice, we take the argument to be the geometric mean 
P = y/PiPA of the two instanton radii. The interaction is attractive for the relative orientation 
cos^ = 1, but vanishes after averaging over all angles 6. The short distance behavior depends 
on the ansatz chosen. In the sum ansatz, there is a substantial repulsive core at distances 
R < \f^p 0, but the amount of repulsion at short distances becomes significantly weaker 
using the more refined trial functions. We will discuss this question in more detail in the 
next section. A parameterization of the interaction in the ratio and streamline ansaetze is 
given in appendix A. In fig. la, we show the ratio and streamline interaction for the most 
attractive and repulsive orientations. One clearly observes that the interactions are similar 
at large distance, but differ significantly at short distance. In particular, the streamline 
interaction has no repulsion at all for the most attractive orientation. The interaction 
smoothly approaches Sint = —25*0 at short distance, corresponding to the annihilation of 
the instanton-antiinstanton pair. Fig. la also shows the effect of a phenomenological core in 
the streamline interaction. The reasoning behind this modification will be discussed in more 
detail in the next section. 
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For all the trial functions except for the simple sum ansatz, the instanton-instanton 
interaction is much weaker than the instanton-antiinstanton one. In fact, in the streamline 
ansatz, the instanton-instanton interaction vanishes. This is a reflection of the fact that 
there is an exact two-instanton solution (with S = 2So) for arbitrary values of the relative 
coordinates. 

The fermionic determinant induces a very nonlocal interaction among the instantons. 
Evaluating this determinant exactly in the instanton ensemble still constitutes a formidable 
problem. In practice we factorize the determinant into a low and a high momentum part 



/N++N. \ 

det{D + ruf) = i Y[lMpi \ det{T + mj), (9) 

where the first factor, the high momentum part, is the product of contributions from indi- 
vidual instantons calculated in gaussian approximation, whereas the low momentum part 
associated with the fermionic zero modes of individual instantons is calculated exactly. 

This means that the instanton induced 't Hooft interaction between quarks P,p7| is 
included to all orders. The low momentum part of the spectrum of the Dirac operator is 
also of special significance in connection with the structure of chiral symmetry breaking. 
TjA is the x matrix of overlap matrix elements 

TiA = j d^x (i)Aoix - ZA)iDx(pio{x - zi), (10) 

where (pi^Ao are the fermionic zero mode wave functions of the instanton and antiinstanton. 
Due to the chirality of the zero modes, the fermionic overlap matrix elements between 
instantons with the same topological charge vanishes. In the following, we will only consider 
quadratic matrices with = A^_. Ensembles with A^_,_ 7^ correspond to systems with 
net topological charge. In the thermodynamic limit, the distribution of winding numbers is 
sharply peaked around zero and imposing an additional condition corresponding to a 
finite winding number is not expected to affect our results. Fluctuations of the topological 



charge can be studied by considering appropriately chosen subsytems |21 



The general structure of the overlap matrix elements is given by Tja = {u ■ R)f{R). This 
means that, like the gauge field induced interaction, the fermionic overlap is maximal when 
the relative instanton-antiinstanton orientation is given by cos^ = 1. Like the gauge field 
induced interaction between instantons, the fermionic overlap matrix elements depend on 
the ansatz for the two-instanton gauge potential. In this case, however, the dependence on 
the trial function is much less pronounced. For the sum ansatz, one can use the equations 
of motion and replace the covariant derivative in (|10]) by an ordinary one. The result can 
be parameterized by 

1 4 

TiA = i{u-R) (11) 

PiPa (2.0 + R^/piPaY 

which is exact at large distances. The streamline ansatz gives the same large distance 
behavior, but somewhat different results at small and intermediate separations. We give a 
parameterization of the streamline matrix elements in appendix B. 



III. THE FREE ENERGY OF THE INSTANTON ENSEMBLE 

In this section, we describe a method to evaluate the partition function of the instanton 
liquid. Using this method, we can calculate the free energy numerically as a function of the 
density of instantons, and determine the equilibrium density from the condition that the 
free energy is minimal. 

The problem in determining the free energy is connected with the fact that the compli- 
cated statistical mechanics associated with the partition function (H) can in general only be 
dealt with by performing Monte Carlo simulations P, p!0| , |39| . These simulations are ideally 



suited for the calculation of various expectation values, but do not give a direct determi- 
nation of the partition function, which provides the overall normalization. Previous Monte 
Carlo calculations have therefore been restricted to simulations of the ensemble at a fixed 
density of instantons, which was determined from phenomenological considerations (typi- 
cally 1 fm~^). Here, we want to go beyond this approximation and minimize the free energy. 
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A method to calculate the partition sum, which is well known in statistical mechanics, quan- 
tum mechanics (see, e.g. and lattice gauge theory [^, is "adiabatic switching". For 
this purpose one writes the effective action as 

Seff = So + aSi, (12) 

which interpolates between a solvable action 5*0 and the full action So + Si. If the partition 
function for the system governed by the action 5*0 is known, the full partition function can 
be determined from 

logZ(a = l) = logZ(a = 0) - da' {0\Si\0)a', (13) 

Jo 

where the expectation value (0|.|0)a depends on the coupling constant a. In our case, the 
effective action is given by 

N++N- 

Seff = - E ^og{d{pi)) + Sint + tTlog{D + mf). (14) 

i=l 

The obvious choice for decomposing the effective action of the instanton liquid would be 
to identify the logarithm of the single instanton distribution with the free action, Sq = 
J2i^og{d{pi)) . This procedure, however, does not work since the instanton distribution 
behaves like d{p) ~ p^^^^\ so that the p integration in the free partition function would not 
be convergent. This is the famous infrared problem which plagues the dilute instanton gas 
approximation [Q. As explained in more detail in the next section, the instanton liquid is 
stabilized by the the repulsive core in the gauge field interaction once the full interaction is 
taken into account. We therefore consider the following decomposition 

N++N- / 2 \ 

Seff= E [~log{d{pi)) + {l-a)u-^) +a{Sint + tT\og{b + mf)) , (15) 

where i' = (6 — 4)/2 and is the average size squared of the instantons with the full 
interaction included. The term proportional to (1 — a) serves to regularize the p integration 
for a = 0. It disappears for a = 1, where the original action is recovered. The specific 
form of this term is irrelevant, our choice here is motivated by the fact that (at least for 
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the one-loop measure d{p)) Seff{a=0) yields a single instanton distribution with the correct 
average size p^. This means that the single instanton distribution generated by S'e//(a=0) 
is a variational ansatz for the full one-instanton distribution. 

The partition function corresponding to the variational single instanton distribution is 
given by 

where /zq is the normalization of the one-body distribution. The p integration in /iq is 



log Z = log(Zo) + N £ da' (0|z/4 - ^ {S^nt + tr log(D + m/)) |0)„., (17) 



regularized by the second term in (|l^). The full partition function obtained from integrating 
over the coupling a is 

where N = + A^_. The free energy density is finally given by F = —1/V- log Z where V 
is the four- volume of the system. The pressure and the energy density are related to F by 

P=-F, ' = T%-P- (18) 

At zero temperature we have e = —p = F and the free energy determines the shift of the 
QCD ground state relative to the perturbative vacuum. Such a shift is certainly present in 
our case, since tunneling lowers the ground state energy. 

The full partition function can be compared to the variational ansatz introduced in |M] 



and employed in many works on the subject |^ , |23|J26| , |24|J35| , |31|| . For simplicity we restrict 
the discussion at this point to pure gauge theory, i.e. neglect the fermionic determinant. 
Since the variational ansatz ignores any correlations between instantons, only the color and 
spatial average of the interaction enters 

r N 

J d^RdUSintiR, U, pi, p2) = ^^'j;^t^pIpI (19) 

where Sint{R,U, pi, P2) is the interaction of two instantons with radii pi^2, separation R 
and relative orientation U. In the sum ansatz, both the II and I A interaction give the 
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same average repulsion = ^vr^ |^. In the ratio and streamline ansatz, this repulsion 
is considerably weaker. In the streamline ansatz, only the I A interaction is repulsive with 
= 4.772 H. 

If the variational single instanton distribution (i(p) exp(— z/p^/p^) is close to the true 
distribution for a = 1 we can calculate the expectation value in (0) using the variational 
one. One finds {Sint) — {Sint)a=o = iVz//2 and the resulting estimate for the partition 
function is 

^ (\/po)^exp(-^), (20) 



((iV/2)!)2^^"^ 2 

which agrees with the result derived in [0. Varying F = —1/V ■ \ogZ with respect to the 
density one finds the expected result NjV = 2fiQ. 



Numerical results for different interactions were compared in P2| , |35| . For the sum ansatz, 
the variational method gives N/V = 0.18A^ with p = 0.48A~^, whereas the streamline ansatz 
gives N/V = 0.54A^ and p = 0.69A~^. The dimensionless diluteness / = 0.57r^p^(A^/l^) 
(the fraction of space-time occupied by instantons) of these ensembles is / = 0.05 and 
0.60 respectively, to be compared with / = 0.06 in the original instanton liquid model. The 
variational method was extended to light quarks in [^,^,^,^ . We will study this problem 
in detail in the next section. 

If correlations among instantons are important, the variational method is not expected 
to provide an accurate estimate for the partition function and other observables. Since 
the main source of correlations in the instanton liquid are dynamical quarks, this issue 
is particularly important for real QCD with two light and one intermediate mass flavor. 
Also, as argued in the introduction, we expect chiral symmetry to be restored due to the 
formation of instanton-antiinstanton molecules. This feature is certainly not captured by 
the variational model (at least not in its simplest form), and we will therefore study the full 
partition function numerically using the method introduced in this section. 
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IV. THE INSTANTON ENSEMBLE AT ZERO TEMPERATURE 



In this section we want to present numerical results obtained from simulations of the 
instanton liquid at zero temperature. While the streamline ansatz provides in principle the 
"best" classical interaction, its derivation relies heavily on conformal symmetry and we do 
not know how to extend it to finite temperature. In section 6, we will therefore give a brief 
discussion of the ratio ansatz ensemble at zero temperature. 

A general problem with calculations in the interacting instanton model is the treatment of 
very close instanton-antiinstanton pairs. On the one hand, these configurations do not con- 
tribute significantly to physical observables such as the quark condensate, hadronic masses or 
the topological susceptibility. Very close instanton-antiinstanton configurations correspond 
to perturbative fluctuations and should not be taken into account as a non-perturbative 
effect. On the other hand, the partition function of the instanton liquid introduced in the 
last section treats even very close pairs as two independent pseudoparticles. Especially in 
the streamline ansatz, which provides very little short range repulsion, this means that close 
pairs can contribute significantly to the free energy of the system. Ideally, one would need 
a consistent determination of the space of collective variables for very close pairs and a sub- 
traction procedure for purely perturbative fluctuations. Unfortunately, such a method is still 
missing. There is an interesting suggestion, defining the instanton interaction (via the opti- 
cal theorem) by the cross section for isotropic multi-gluon production [0]. In this case, the 
existence of a short range core in the instanton interaction is related to the question whether 
the multi-gluon production cross section at high energy can reach the unitarity bound. If, 
as it has been suggested in |^ and other works on the subject, the cross section only grows 
until it reaches a square root suppression and then turns, the instanton-antiinstanton inter- 
action would decrease until it reaches a minimum at which Sint = —Sq, and then show a 
short distance repulsive core. 

In practice we have decided to deal with the problem of very close pairs in the streamline 
ansatz by introducing a purely phenomenological short range repulsive core 
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in both the II and lA interactions. Here A is the conformal parameter that determines 
the functional form of the streamline interaction P5| and the parameter A controls the 



strength of the core. This parameter essentially governs the dimensionless diluteness / = 
0.57r^p^(A^/l^) of the ensemble. The second parameter of the instanton liquid is the scale 
Aqcd in the partition function, which fixes the absolute units. Although Aqcd is known 
in principle from applications of perturbative QCD, the accuracy of these determinations is 
not very high and quantities like the instanton density (N/V) ~ Aq^^ are very sensitive 
to the precise numerical value of the scale parameter. As in lattice calculations, one may 
therefore fix the scale using the value of some observable, for example a hadronic mass, 
as the basic unit. Which quantity to hold fixed while comparing different theories, e.g. 
quenched and unquenched calculations, is pure convention. We have decided to proceed in 
a very simple way and fix our units such that N/V = 1 fm~^ in all cases. This means that 
in this work Ifm is, by definition, the average distance between instantons (at T = 0). The 
corresponding instanton density agrees well with the lattice measurements mentioned above 
and corresponds to the canonical value of the gluon condensate, (^G^) = (350 MeV)^. 

We have studied the instanton ensemble for various values of the core parameter A. The 
results are summarized in table 1. By increasing the core one obviously makes the ensemble 
more dilute. However, as the interaction becomes more repulsive, the density in units of 
the scale parameter drops and one has to increase the value of Aqcd in order to keep the 
physical density fixed. In practice, we have chosen A = 128, so that Aqcd remains below 300 
MeV, which is roughly the upper limit of the experimental uncertainty. As a consequence, 
the ensemble is not as dilute as suggested by phenomenology (n = 1 fm~^ and p = 0.33 fm), 
but comparable to the lattice result n = (1.4-1.6) fm~^ and p = 0.35 fm [|r^. Two more 
important parameters of the ensemble are the average instanton action S = {^'k'^)/ and the 
average interaction Smt/N per instanton. For A = 128 we find 5* ^ 6.4 and Smt/N ~ 1.0, 
showing that the system is still semiclassical and that interactions among instantons are 
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important, but not dominant. 

Detailed results of our simulations for A = 128 are shown in fig. 2-5. The partition 
function at each instanton density was determined by generating 5000 configurations with 32 
instantons at 10 different coupling constants a. The variational ansatz (|16|) for the partition 
function was determined from 600 initial sweeps with the full interaction {a = 1). The 
integral ([17|) was determined by gradually lowering the coupling to a = and then raising 
it back to a = 1 (hysteresis method). The difference in the result between the up and down 
sweeps provides an estimate of the error in the integral due to incomplete equilibration. The 
average between the up and down sweeps usually provides a good estimate for the correct 
result, even if equilibration is slow (as will be the case close to the phase transition). 

Fig. 2a shows the free energy versus the instanton density (in units of A^) for the pure 
gauge theory (without fermions). At small density the free energy is roughly proportional 
to the density, but at larger densities repulsive interactions become important, leading to 
a well-defined minimum. We also show the average action per instanton as a function of 
density. The average action controls the probability exp(— S*) to find an instanton, but has 
no minimum in the range of densities studied. This shows that the minimum in the free 
energy is a compromise between maximum entropy and minimum action. 

Fixing the units such that N/V = 1 fm~^, we have A = 270 MeV and the vacuum energy 
density generated by instantons is e = — 526 MeV/fm'^. This important quantity is related 
to the gluon condensate by the trace anomaly 

where b = yA^'c in the pure gauge case. For a sufficiently dilute system of instanton^, 
the gluon condensate is simply proportional to the instanton density, e = —b/A{N/V) = 
-565MeV/fml The good agreement of this number with the energy density determined 



^As mentioned above, the average interaction is only a 15% correction to the single instanton 
action. 
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above shows that our calculation is consistent with the trace anomaly. For the variational 
calculation, this is also the case as indeed one would expect for any calculation that does 
not introduce any dimensionful scale except for A. Also note that not only the depth of 
the free energy, but also its curvature (the instanton compressibility) is fixed from a general 
theorem. From standard thermodynamical relations, the compressibility is related to the 
mean square fluctuation of the particle number. This quantity is fixed from a low energy 
theorem derived by Novikov et al. |^ 



d^x < g'G\x)g^G\0) > = < g'G' > . (23) 

Similar to the trace anomaly (0), this relation follows from the renormalization properties of 
QCD. Saturating (^31) with a dilute system of instantons, one has ((AA^)^) = |A^, where AA^ 
is the mean fluctuation of the number of instantons in a volume V. For the compressibility 
of the instanton liquid this implies 

d{N/vy ~ 4 \ v) ■ 

We have determined the compressibility from our data and find 3.2{N/V)^^, to be compared 
with 2.75(A^/l^)~^ from eq.(^). We have also studied density fiuctuations in the interacting 
instanton liquid pT|JTB[] and found good agreement with the low energy theorem (see also 



(24) 



MM)- 

Next we study the effects of light quarks on the instanton ensemble. A problem that is 
well known in lattice simulations is the fact that finite volume effects are more severe in the 
presence of light fermions. These "mesoscopic" effects become important if the mass of the 
Goldstone bosons ~ ^/niA is comparable to the inverse box size 1/L. In order to avoid 
these problems one is often forced to work with light quark masses that are larger than their 
physical values. The calculations reported below were performed with light quark masses 
rriu = rrid = 0.1 A, while the strange quark mass is = 0.7A. If one makes the quark mass 
smaller then chiral symmetry breaking will eventually disappear due to finite size effects. 

The free energy as a function of the instanton density in full QCD is shown in fig. 3a. 
We find that the free energy looks very similar to the pure gauge case, but the minimum 
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is shifted to a smaller density N/V = 0.174A"'. With our convention N /V = lfm~^ as 
above, the scale is given by A = 310 MeV. This gives the instanton-induced vacuum energy 
e = — 280 MeV/fm'^, which is smaller as compared to the pure gauge case. In the presence 
of quarks, the trace anomaly gives[| 

' = "12^^ ^ > +E ^"^/ < Wf > ■ (25) 

The main difference as compared to the pure gauge case is that the Gell-Mann-Low coefficient 
is changed from 6=11 to 6 = 9 (for Nf = 3). Estimating the gluon condensate from the 
instanton density as above, we get e = — 488MeV/fm^. This number does not agree very 
well with the value extracted from the free energy. One possible reason is that interactions 
are more important in the full ensemble so that the gluon condensate is not proportional 
to the instanton density. Another problem might be connected with the value of the quark 
mass. As mentioned above, if one enters the regime of mesoscopic QCD, one may encounter 
additional scale breaking not described by the trace anomaly. We indeed find that the scale 
anomaly can be saturated with the instanton density for larger quark masses. 

At the minimum of the free energy we can also determine the quark condensate (see 
fig.3c). In full QCD we find (qq) = — (216MeV)^, which is in good agreement with the phe- 
nomenological value. This value can also be compared to the quark condensate in quenched 
QCD (calculated in the pure gauge configurations, see fig. 2c) {qq) = — (251MeV)^, showing 
that light quarks suppress the quark condensate. 

Additional information about the ensemble is provided by the distribution of the instan- 
ton sizes, Dirac eigenvalues and fermionic overlap matrix elements. In figures 4 and 5 we 
present these distributions for the two cases discussed above. The instanton size distribu- 
tion shows the perturbative p^~^ behavior at small sizes, has a maximum (at p = 0.50A~^ 

^Despite the factor 1/4 in front of the quark condensate contribution, this relation is consistent 
with the general result {qfqf) = log Z. The remaining mass dependence of the vacuum energy 
density comes from the ghion condensate, g^{9^G^) = ~^^f-{<if<lf) ^ which is valid for light quarks. 
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in quenched and p = 0.65A~^ in full QCD) and falls off for large sizes. The average sizes 
are p = 0.43 fm in quenched and p = 0.42 fm in full QCD. The diluteness parameter in the 
full ensemble is / = 0.14, to be compared with / = 0.17 in the quenched calculation. 

The distribution of eigenvalues of the Dirac operator, Dip\ = Xip\, contains a lot of useful 
information about the spectrum of the theory. We will discuss some of these questions in 
more detail below. Here we only mention that the Banks-Casher relation (qq) = — p(A = 
0)/7r connects the density of eigenvalues of the Dirac operator near zero with the chiral 
condensate. From figs. 4b and 5b one clearly observes that the instanton liquid leads to a 
nonvanishing density of eigenvalues near A = 0. The results also show that the presence 
of light quarks suppresses the number of small eigenvalues. While the spectrum is peaked 
towards small eigenvalues in the purge gauge case, it is essentially flat in full QCD. This is 
consistent with the prediction from chiral perturbation theory |^ 



which is valid for Nf > 2 (speculations about the Dirac spectrum in quenched QCD can be 
found in ^^). The second term is zero0 for Nf = 2 and leads to a singular behavior for 
Nf > 2. No such dip is seen in our spectrum, showing that the result is closer to the Nf = 2 
case. More detailed studies of the spectrum of the Dirac operator in the instanton liquid for 



various number of colors and flavors were reported in [50 



In fig. 3c and 4c we also show the distribution of the largest fermionic overlap matrix 
elements Tja_ for each instanton. For each instanton we select the antiinstanton it has the 
largest overlap with and plot the resulting distribution of matrix elements. This distribution 
is a measure of the strength of correlations among the instantons. For a completely random 
system at the same density as the simulated quenched and fully interacting ensembles, 
the average overlap would be Tja = 0.22A and Tja = 0.21A, respectively. Instead, the 



^This result is connected with the fact that for Nf > 2 there is a Goldstone boson cut appearing 
in the scalar isovector (5 meson) correlator, while there is no (5 — > vrvr decay allowed for two flavors. 
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measured distribution have an average 0.23A and 0.33A, showing that correlations are not 
very important in the quenched ensemble, but play some role in the full ensemble. A 
more physical measure of the importance of correlations among instantons is given by the 
dependence of hadronic correlation functions on the different ensembles. We will study this 
question in some detail in a forthcoming publication. 



V. THE PHASE DIAGRAM OF THE INSTANTON LIQUID 

After fixing the the overall scale at zero temperature it is now straightforward to ex- 
tend our calculation to finite temperature. In a euclidean field theory, finite temperature 
only enters through the boundary conditions obeyed by the fields. The gauge fields have 
to be periodic in the euclidean time direction with the period given by the inverse temper- 
ature, while fermions are subject to antiperiodic boundary conditions. The corresponding 
periodic instanton and fermion zero mode profiles can be constructed using 't Hooft's multi- 
instanton solution |51[]. These profiles may then be used to study the instanton interaction 
and fermionic overlap matrix elements. For the ratio ansatz, such a study was performed in 



5^ . A parameterization of this interaction is given in appendix C and D. 

The most important qualitative feature of the instanton interaction at finite temperature 
is the form of the fermionic overlap matrix elements. The determinant for one instanton- 
antiinstanton pair separated by r in euclidean time and r in the spatial direction is roughly 
proportional to 



det(D) 



sm(7ri T) 



(27) 



cosh(7rTr) 

The form of this interaction is a simple consequence of the periodicity in the time direction 
and the fact that fermion propagators are screened in the spatial direction by the low- 
est Matsubara frequency vrT. The interaction clearly favors instanton-antiinstanton pairs 
(molecules) that are aligned along the time direction with a separation r = 1/(2T). II- 
genfritz and Shuryak proposed that this feature of the interaction leads to a phase 
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transition in the instanton liquid. In this phase transition, the system goes from a random 
hquid to an ordered phase of instanton-antiinstanton molecules. The transition was studied 



in a schematic model in [^,^]. In the present work we want to test whether a consistent 
treatment of the full partition function of the instanton liquid does indeed lead to the ex- 
pected phase transition and perform a quantitative study of bulk properties of the system 
near the transition. 

Before we consider the QCD case (with two light and one intermediate mass flavor) in 
detail, it is instructive to study the phase diagram of the instanton liquid for a wider range 
of theories. In order to cover many different parameters, we have restricted ourselves to an 
exploratory study, in which we do not minimize the free energy of the system, but consider 
the instanton ensemble at fixed density. Using the same convention employed in the last 
section, this density was chosen to be Ifm""^. In our experience, if the system undergoes a 
phase transition, it will do so even if the density is not determined self-consistently, although 
at a somewhat different temperature. Nevertheless, before this point is checked in detail, 
the results presented in this section should be considered qualitative. 

We first consider the situation for only one quark flavor. In this case the only chiral 
symmetry is the anomalous axial f/^(l), so there are no Goldstone pions, but only a flavor 
singlet (the analog of the t]'). Therefore, one would not expect a chiral restoration phase 
transition. On the other hand, a (first order) transition may exist even without a symmetry, 
and in principle, the general argument concerning the formation of I A molecules also applies 
to the case Nf = 1. Clearly, fermion induced correlations become stronger as the number of 
flavors increases, so whether they are strong enough for Nf = 1 to induce a phase transition 
has to be studied by performing simulations. 

We have studied the Nj = 1 system at an instanton density N/V = lfm~^ for tempera- 
tures up T ~ 300 MeV. We have observed no phase transition in this range, the condensate 
(qq) decreases smoothly and is non-zero even at the highest temperature studied. Of course, 
even in pure SU(3) Yang Mills theory (without fermions) lattice simulations find a strong 
first order deconfinement phase transition at a temperature T ~ 250 MeV (where the scale 
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is set by the string tension or quenched p meson mass). The interacting instanton model, 
however, does not give confinement, so, naturally, there is no deconfinement transition. 
The case of two light flavors, Nj = 2, is quite special. According to standard universality 



arguments |]53| , |5^ , the chiral phase transition is expected to be of second order, with the same 
critical behavior as the 0(4) Heisenberg magnet in d = 3 dimensions 0. We have simulated 
Nf = 2 ensembles for various temperatures and values of the quark masses. For sufficiently 
small quark masses m < 0.15A we find large fiuctuations in the condensate at T ~ 150MeV. 
To illustrate this fact, we show a typical time history of the quark condensate in fig. 6a. Every 
configuration corresponds to one complete sweep through the instanton ensemble, with one 
Metropolis hit performed on every collective variable. It is clear that the transition is either 
second order or first order with a very small barrier between the two phases. In order to 
distinguish these two possibilities one would have to perform a finite size scaling analysis, 
which goes beyond the scope of this paper. We will see later that for more fiavors there is 
a barrier between the two phases. In fig.6b we show the Dirac spectrum for T = 150 MeV 
and m = 5 MeV. The eigenvalue density extrapolates to zero at A = 0, but the slope is very 
steep and there are no signs of a gap in the spectrum. 

We will treat the three fiavor case in detail in the next section, and now proceed to the 
case of several fiavors of light quarks. Let us first mention several theoretical arguments 
suggesting qualitative changes in the vacuum structure as Nf increases beyond some critical 
value Np\ 

• A simple phenomenological observation is the fact that for Nf > 12 the number of 
"pions" Nt, = Nf — 1 exceeds the number of quarks and antiquarks Nq = ANcNf. This 



^Recently, Kocic and Kogut [^] have challenged these arguments and suggested that the transition 
is indeed second order, but with mean field exponents. Lattice results (see and the review |]57| ) 
seem to verify that the transition is second order, but are not yet sufficiently accurate to resolve 
the issue about the critical exponents. 
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prevents the usual matching of the pion gas at low temperature to the high-T quark- 
gluon plasma phase with a positive bag constant. It seems plausible that something 
should happen before this point is reached, so A^^*"** < 12. 

As noticed by a number of authors (see, e.g. the second coefficient V of the beta 

function changes sign for smaller N j than the first one, h. In this window one has the 
possibility of asymptotic freedom coexisting with an infrared fixed point. 



It was recognized in that for large N ^ the instanton liquid is dominated by 

instanton molecules, and that their density is UV divergent (at small radii) if 26— 5 < 0, 
corresponding to N j > 12 (this estimate is slightly modified if interactions are taken 
into account ^^). This divergence, however, is a short distance effect and should not 
affect physical observables. 



• Lattice simulations of QCD for large Nf were performed by Iwasaki et al. |£0[- They 
found that for Nc = 2 the critical number of fiavors for which chiral symmetry is 
restored is N^"^ = 3, while for A^^ = 3 it is Nj^^^ = 7. A strange "bulk transition" was 
also observed by the Columbia group for Nf = 8 |6T| . 

• Finally, significant progress has been made in understanding A^ = 1 supersymmetric 
generalizations of QCD ||62|. One result that is of particular interest in the present 
context is the exact determination of the critical number of fiavors NJ"^ = A'c + 1, 
where chiral symmetry is broken for Nf < Nf"^, but not for Nf > Nf"^. 

Let us now report on our results concerning the vacuum structure of the instanton liquid 
for many fiavors. A useful tool to analyze the structure of chiral symmetry breaking that 
was suggested by the Columbia group is the "valence quark mass" dependence of the 



quark condensate, defined by 

2m^ 



<qq{m,)> = - j dX p{X,m)^^-^ . (28) 



Here, p(A, m) is the eigenvalue density of the Dirac operator calculated from configurations 
generated with dynamical mass m. The quark condensate in the chiral limit is defined 
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by going to the thermodynamic hmit and then taking m = rriy — > 0. In a finite system, 
however, chiral symmetry is never broken at nonzero quark mass and one has to perform a 
detailed scahng analysis. This analysis amounts to performing many simulations at various 
values of the dynamical mass, which is a very time consuming procedure. Instead, one 
might get a good indication of the full behavior by studying the valence mass dependence 
of the condensate. At very small m^, the condensate is proportional to m^, which means 
that one is sensitive to induced rather than spontaneous chiral symmetry breaking. At very 
large mass m^, the condensate is proportional to m~^. Both of these limits are of course 
unphysical, and spontaneous symmetry breaking in the continuum limit is indicated by the 
appearance of a plateau at intermediate values of the valence mass where (qq) depends only 
weakly on m„ . 

In fig. 7 we show a number of calculations of the valence mass dependence of the conden- 
sate. The results in fig. 7a were obtained for 64 instantons in a cubic box V = (2.828 fm)^ 
(corresponding to a density N/V = lfm~^ and a fairly low temperature, T = 71 MeV) 
with the dynamical quark mass m = 20 MeV and different number of flavors Nf = 1, . . . , 8. 
Regions where the condensate depends only weakly on rriy are clearly seen for Nf = 1,2,3. 
These plateaus are clearly absent for Nf = 5 or larger: we therefore conclude that the in- 
stanton liquid model has a chirally symmetric ground state in these cases. A more detailed 
study of the conflgurations shows that all instantons are bound into molecules. The border- 
line case appears to be Nf = 4: if a condensate is present, it is signiflcantly smaller than in 
QCD and due to a relatively small random component of the vacuum. 

We have also performed simulations for two colors, and obtained a similar picture, but 
with the critical number of flavors smaller by one. For A^^c = 2 no clear signal of chiral sym- 
metry breaking was observed for Nf > 3. This suggests the empirical relation Nj"^ = Nc + C 
with 1 < C < 2, an amusing coincidence with Seiberg's results for A^ = 1 supersymmetric 
QCD 0. 

The next question we would like to address is the dependence of these results on the 
dynamical quark mass m. Increasing the dynamical mass, one decreases the influence of 
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quark determinant. A large quark mass works against the correlations induced by a large 
number of flavors. This means that for a sufficiently large quark mass one can flnd "spon- 
taneous" symmetry breaking even if there is no symmetry breaking in the chiral limit. This 
is demonstrated in Fig.7b, where we show a series of calculations for Nf — 5 with the dy- 
namical quark mass m ranging from 20 to 100 MeV. In this case a plateau in the valence 
mass dependence of the condensate reappears at rather small critical quark mass of about 
35 MeV. 

Repeating these studies for different numbers of flavors and varying the mass and tem- 
perature, one can map the phase diagram of the instanton liquid in the T — m (temperature- 
mass) plane. In flg.8 we have summarized our results for Nf — 2,3 and 5. Those cases in 
which we found a signal for chiral symmetry breakdown are marked by the open squares, 
and those corresponding to the restored phase by solid squares. The (approximate) location 
of the discontinuity line between the two phases is marked by the stars connected by dashed 
lines. For two flavors we do not flnd such a line of discontinuities. We have performed 
a number of runs in the vicinity of the phase transition, all showing large fluctuations of 
the condensate. For larger Nf we clearly see a transition line with a discontinuity of the 
condensate. When the number of flavors increases, one end of this line moves to the left (the 
critical temperature Tc decreases with Nf), crossing zero somewhere around Nf = 4, but 
before Nf — 5. In these cases, the ground state exhibits spontaneous symmetry breaking 
only if the quark mass exceed some critical value. For technical reasons we have not tried 
to follow the discontinuity lines much beyond T > 200 MeV. In this regime, the condensate 
becomes more and more dominated by very few extremely small eigenvalues, so that numer- 
ical accuracy starts to become a concern. In addition to that, one does not expect to be 
able to describe possible transitions at very large temperature in the instanton model, since 
these will presumably be dominated by the pure gauge deconflnement transition. 
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VI. THE INSTANTON ENSEMBLE IN QCD AT FINITE TEMPERATURE 



In this section we want to discuss in detail the physically relevant case of two light and 
one intermediate mass flavor. In particular, we will perform self consistent simulations where 
the correct instanton density is determined from minimizing the free energy. The instan- 
ton interaction and fermionic matrix elements that enter these calculations have already 
been discussed in the last section. The semiclassical calculation that leads to the instanton 
distribution (|^) has also been generalized to finite temperature giving 

d{p, T) = dip, T = 0) exp(-l(2iV, + Nf) [ixpTf - Biz)) (29) 

with z = TipT. This correction factor leads to an exponential suppression of large instantons 
at high temperature. Its origin is mainly the scattering of thermal gluons on the instanton. 
This phenomenon is related to the Debye mass, which in fact has the same dependence on 
Nc and Nf. The applicability of ( ^91) is controlled by two separate conditions, p <^ hr^ 
to ensure the semiclassical treatment, and T ^ A to justify a perturbative treatment of 



the heat bath. It was therefore argued in p8[ that one should not use the perturbative 



suppression factor ( pOf ) at temperatures below the phase transition. This suggestion was 
indeed verified by lattice simulations for pure gauge theory [^], which found a very week 



temperature dependence of the instanton density below T^, and an exponential suppression 
of instantons consistent with (EUf) above T^.. In practice we have determined the phase 



transition temperature without the suppression factor (EUp . The final simulations were then 



performed with an additional temperature dependent factor (1 — tanh(2^^))/2 inserted in 
the exponential appearing in eq. (^Of) . From our simulations, we have determined — 0.75A. 
We have taken AT ~ 0.2A. This number characterizes the temperature range over which 
the Debye mass reaches its perturbative value and is roughly consistent with the result 
of lattice simulations [Q. In our simulations, the value of AT determines how fast the 
instanton density drops above Tc. On the other hand, we have checked that the temperature 
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dependence of the quark condensate is almost independent of this parameter. 

With all the ingredients of the partition function fixed, we can proceed to our simulations. 
The free energy is determined as described in sections 3 and 4. There are no difficulties 
in principle associated with applying the coupling constant integration method at finite 
temperature. In practice, however, some care is required since in case the system undergoes 
a phase transition, the ensemble at the full coupling a = 1 is in a different phase as compared 
to the random distribution at a = 0. This implies that the transition will occur at some value 
of a during the coupling constant integration. The hysteresis method described in section 
4 is particularly well suited to handle a situation like this. Physically, this simply means 
that the variational ansatz ([TBD , while it is still used as a reference distribution, ceases to 
be a good approximation to the fully interacting instanton distribution. In the whole range 
of temperatures studied, we find that the behavior of the free energy as a function of the 
instanton density looks qualitatively similar to the zero temperature result shown in figs. 2 
and 3. In particular, the free energy is negative at small density and has one fairly well 
defined minimum. This means that the phase transition is either smooth or a first order 
transition with a small jump in the instanton density, which we cannot resolve due to finite 
size effects. 

The resulting instanton density, free energy and quark condensate are shown in fig. 9. In 
the ratio ansatz the instanton density at zero temperature is given by N jV = 0.69A^. Taking 
the density to be 1 fm~^ at T = fixes the scale parameter A = 222 MeV and determines the 
absolute units. The temperature dependence of the instanton density is shown in fig. 9a. It 
shows a slight increase at small temperature^, starts to drop around 115 MeV and becomes 
very small for T > 175 MeV. The free energy closely follows the behavior of the instanton 
density. This means that the instanton-induced pressure first increases slightly, but then 



The zero temperature point corresponds to a simulation with zero temperature matrix elements 
in a periodic box V = (2.828A^^)^, so it is not a true zero temperature calculation. 
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drops and eventually vanishes at high temperature. This behavior is expected for a system 
of instantons, but in a complete theory with perturbative effects, the pressure should always 
increase as a function of the temperature. In QCD, the pressure at high temperature is of 
course provided by the black body contributions from quarks and gluons. We will come 
back to this issue at the end of the section. 

The temperature dependence of the quark condensate is shown in fig. 9c. At temperatures 
below T — 100 MeV it is practically temperature independent. It then starts to drop fast 
and becomes very small around the critical temperature T ~ 140 MeV. Note that at this 
point the instanton density is N/V = 0.6 fm~"^, slightly more than half the zero temperature 
value. This means that the phase transition is indeed caused by a transition within the 
instanton liquid, not by the disappearance of instantons. This is illustrated by fig. 10, which 
shows projections of the instanton liquid at temperatures T = 74, 123, 158 MeV, below, near 
and above the chiral phase transition. The figures are projections of a 4-dimensional cube 
V — (3.00A~^)^ • 1/T into the 3-4 plane. The positions of instantons and antnnstantons are 
denoted by +/— symbols. The lines connecting them indicate the strength of the fermionic 
overlap matrix elements. Below the phase transition, there is no clear pattern. Instantons 
are unpaired, part of molecules or larger clusters. As the phase transition progresses, one 
clearly observes the formation of polarized instanton-antiinstanton molecules. 

More details of the phase transition can be inferred from the graphs in figs. 11, 12 and 13. 
Fig. 11 shows trajectories of the quark condensate for the temperatures mentioned above. 
Every configuration corresponds to one complete sweep through the instanton liquid, with 
one Metropolis hit on the collective variables of all instantons. As discussed in more detail in 
the next section, the fiuctuations in the quark condensate determine the scalar susceptibility. 
Below the phase transition, these fluctuations are controlled by the effective a meson mass. 
Near the phase transition, we find evidence for a weak first order transition with the system 
frequently tunneling between the two phases. In order to clearly distinguish this scenario 
from a smooth crossover we would have to run with significantly larger volumes. Fig. 12 
shows the spectrum of the Dirac operator. Below the phase transition it has the familiar 
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flat shape near the origin and extrapolates to a nonzero density of eigenvalues at A = 0. 
Near the phase transition the eigenvalue density appears to extrapolate to as A —> 0, but 
there is a spike in the eigenvalue density at A = 0. This spike contains the contribution from 
unpaired instantons. Near the phase transition, some of these instantons come from the 
breakup of molecules as the systems tunnels from the high to the low temperature phase. In 
addition to that, even at very high temperature one expects a finite concentration 0{m^f) 
of random instantons. We will further comment on this issue in the next section. 

If there is evidence that the transition is weakly first order, it is of interest to a least 
have an upper bound on the latent heat associated with the transition. We have already 
mentioned that we do not directly observe a latent heat from a discontinuity in the instanton 
density. Another way to calculate the latent heat is using the Clausius-Clapeyron relation 
generalized to QCD 

disce = Tc \ — -] disc < qq > . (30) 

As Leutwyler already realized in his original paper, this relation is useful in analyzing nu- 
merical simulations even if one does not directly observe a jump in e. We have determined 
the derivative of the critical temperature with respect to the quark mass from simulations 
at several different quark masses (see section 7). We find dTc/drriq ^ 0.9 at nig = O.IA, 
close to Leutwyler's estimate dTc/drrig ^ 1. From the trajectories of the quark condensate 
(fig.ll) one concludes that the discontinuity is at most disc(gg) < 0.4(gg)r=o- These num- 
bers give disce < 60MeV/fm^, which is quite small as compared to the energy density at 
zero temperature. Saturating the trace anomaly (e — 3p)/4 = —j{N/V) with instantons, 
one can also translate this estimate into a bound on the jump in the instanton density, 
disc (N/V) < 0.03 fm~'^, to be compared with the zero temperature value (N/V) = lfm~^. 
This bound shows that we do not expect to see an appreciable jump in the instanton density, 
even if there is a first order transition. 

In fig. 13 we show the distribution of cos^(a) = |'U4p/|Mp, the relative color orientation 
angle between instantons and antiinstantons. For this purpose we have selected for every 
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instanton the antiinstanton it has the largest overlap with and calculated the corresponding 
orientation angle cos^(q;). The resulting distribution was normalized to the SU{3) invariant 
measure. Below the phase transition there is no preferred direction and the distribution is 
flat (the spike near cos^(a) = 1 comes from the fact the measure goes to zero in this region). 
Near the phase transition the distribution is strongly peaked towards cos^(a) = 1, showing 
that instanton antiinstanton pairs are indeed strongly polarized. 

Completing this section, let us add a few general remarks concerning the thermodynamics 
of the instanton liquid. In this paper, we have only determined the free energy associated 
with instantons. Our calculation includes the contributions from the zero mode determinant, 
which describes the excitation of low energy collective modes (the pions). However, it 
does not include non-zero mode contributions from quarks or gluons. In particular, it does 
not include the Stefan-Boltzmann contribution at large temperature, or perturbative 0{as) 
corrections to it. In a more realistic description of the thermodynamics of the chiral phase 
transition, these contributions should certainly be included. (An attempt to do so was made 
in the schematic cocktail model by Ilgenfritz and Shuryak |]^.) In this paper, however, it was 



our intention to study the physics contained in the partition function of the instanton liquid, 
without making further assumptions or inputs. Remarkably, the instanton liquid undergoes 
a phase transition in the expected temperature regime even if no additional ingredients are 
added. 

VII. DIRAC EIGENVALUES AND SUSCEPTIBILITIES 



Beginning with the classical paper by Banks and Casher [Q, many useful relations 



have been pointed out between various observables and the spectrum of the Dirac operator 
Dipx = A^/'a- In quantum field theories the region of small A ^ is the analog of the Fermi 
surface in solid state physics, and understanding the density of states near this point is 
crucial for many properties of the system, both at zero and non-zero temperatures. 

Before we discuss our data in more detail we would like to briefly review some of the 
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results that can be found in the hterature. As T — > Tc the quark condensate vanishes. If 
there is a second order phase transition for Nf — 2, the mass dependence of the condensate 
is governed by the critical index S 

<qq>L-T ^m^'\ (31) 

where T^^ax denotes the pseudocritical temperature, corresponding to the position of the peak 
in the scalar susceptibility (see below). As mentioned above, there is still some controversy 
concerning the values of the critical indices. According to standard universality arguments, 
QCD with Nf = 2 case is analogous to the 0(4) Heisenberg magnet, which has 1/5 = 0.21. 
On the other hand, mean field scahng would give 1/S — 1/3. 
Using the Banks-Casher relation 

/2Tn 
dXp{X,m,T)^^^ (32) 

one can try to convert this relation into a prediction for the spectral density of the Dirac 
operator. The problem is that there are two sources for the mass dependence of the con- 
densate, the explicit mass term in the integral (the "valence mass") and the implicit mass 
dependence of the spectral density (the " sea mass" ) . Combining these two effects, one can 
only conclude that 

p(A,T = T^„,)^Ai/^-V, (33) 

where k, characterizes the explicit mass dependence of the spectral density. Still, one may 
hope that the influence of the sea mass is small and study the behavior of the spectral density 
near zero at the phase transition. In fig.6 we have shown the distribution of eigenvalues at 
Tc for the case of two flavors. The spectrum behaves roughly like p(A) ~ A*^'^. In fig. 12 
we have displayed the analogous result for QCD. Here, the spectrum consists of two parts, 
a smooth part behaving like p(A) ~ A^ and a spike near zero. There are two effects that 
contribute to the spike near zero. One is the fact that near a first order phase transition 
both the broken and the restored phase contribute to the spectrum. The other is that at 

30 



finite mass there will always be a finite 0{m^f) number of unpaired instantons, giving a 
contribution of the form p(A) ~ m^f6{X) to the Dirac spectrum. 

Additional information about the phase transition is provided by mesonic susceptibilities. 
We define these susceptibilities as the integral of the corresponding mesonic correlation 
function 

Xr = J d^x < qix)Tq{x)q{0)Tq{0) >, (34) 

where F is a spin-isospin matrix with the appropriate quantum numbers. They characterize 
the response of the system to slowly varying external perturbations. For example, the scalar- 
isoscalar susceptibility can also be defined as the second derivative of log Z with respect to 
the quark mass. Near a second order phase transition, the susceptibilities associated with 
order parameter fluctuations are expected to diverge in the chiral limit in a universal manner. 
In particular, we have the predictions ||5^ 



= 1/5. (35) 

T=T, 



XlT 

Above the phase transition, chiral symmetry is restored and one has Xo- = Xn as the quark 
mass is taken to zero. Below the phase transition, Xcr/Xn —>■ as m ^ 0, since Xn has a 
l/m singularity in the chiral limit, while Xo- only has a logarithmic singularity. As noted 
by Karsch and Laermann [Q, the second relation in ( p5| ) is very useful in determining the 



exponent 6, since measurements of A(m, T) = Xa/Xn at different masses are expected to 
cross at a unique temperature T = Tc, where the value of 1/5 can be extracted from A(m, Tc) 
(up to corrections from the non-singular part of the free energy). 

The susceptibilities associated with flavored mesons like the pion and 6 meson are directly 
related to the spectrum of the Dirac operator. Inserting the general decomposition of the 
quark propagator in terms of eigenfunctions into the definition (0), one finds 

Xn = 2 J d\p{X,m,T)j^^—^, (36) 



+ m 
— 

(A^ + m 



X5 = '2 I dXp{X,m,T)j^^ —. (37) 
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Comparing (|36D with the Banks-Casher relation (|32D one notes that Xn = — (qq) /fn, a rela- 
tion that can also be obtained by saturating the pion correlator with one-pion intermediate 
states and using PCAC The susceptibility in the scalar-isoscalar {a meson) channel receives 
an additional contribution from disconnected diagrams, which can not be expressed in terms 
of the spectral density 



Xa = X5 



+ 2V (< {qqf >-<qq>^). (38) 



The second term measures the fluctuations of the quark condensate in a finite volume V 
(as the volume is taken to infinity). Note that all the susceptibilities (|36| - p8D suffer from 



ultraviolet divergencies related to the behavior of p(A) at large A. Here we are interested in 
infrared divergencies that occur when the quark mass is taken to zero, related to the spectral 
density at small virtuality A. This means that the critical behavior is unrelated to the 
ultraviolet behavior which has to be regularized in some fashion. In our case, the eigenvalue 
density is calculated in a basis spanned by the zero modes of the individual instantons. 
This means that the high virtuality part of the spectrum that causes the susceptibilities to 
diverge is not included anyway. 

One important aspect of the QCD phase transition is the question whether the U{1)a 
symmetry is restored during the transition If this is the case, one expects (for Nf = 2) 
two additional degrees of freedom to become light during the transition, the rj' and 6. The 
rj' susceptibility is connected with fluctuations of the number of zero modes of the Dirac 
operator and cannot be measured directly in a system with total topological charge zeroQ. 
The 6 susceptibility, however, can easily be studied in our simulations. If U{1)a symmetry 
is restored, we expect xs to peak at the transition and diverge as the quark mass is taken 
to zero. Furthermore, above the transition one would have Xtt = X<5 ^is m ^ 0. This also 
requires that the disconnected part of Xa has to vanish for T > as m — > 0. 



^It was demonstrated in |21| that one can study the topological susceptibility by studying fluctu- 
ations of the topological charge in subvolumes of a large system with total charge zero. 
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We show our results for the mesonic susceptibihties in fig. 14 and 15. The results in fig. 14 
were obtained for three fiavors with masses m„ = rrid = O.lOA = 22 MeV and = 0.70A = 
155 MeV. They correspond to the instanton density and free energy discussed in the last 
section. The susceptibilities were obtained from 1000 configurations separated by 5 sweeps 
through an ensemble of 64 instantons. In our simulations, the physical volume depends on 
the temperature, since the number of instantons is held fixed while their density drops with 
temperature. Near the phase transition, we have V3 = (3.5 fm)^. In fig. 15 we show our 
results for a slightly smaller light quark mass m„ = = 0.07A = 15 MeV. 

Fig. 14 clearly shows a peak in the sigma susceptibility, indicating the chiral phase tran- 
sition. The position of the peak is at T ~ 125 MeV, somewhat too low as compared to the 
estimate Tc = 140 MeV from lattice simulations. The delta susceptibility does not show any 
pronounced enhancement. This is also seen from the lower part of the figure, where we show 
the disconnected part of the scalar susceptibility and the difference ~ Xs- Clearly, the 
peak in Xa comes completely from the disconnected part. Furthermore, Xn — Xs becomes 
very small above the phase transition, but shows no tendency to go to zero in the range of 
temperatures studied. At the smaller mass (see fig. 15) there is a broad enhancement visible 
in Xs and x-k — Xs appears to vanish above the phase transition. On the other hand, if one 
compares the peak heights of Xa for the two different quark masses, one finds Xcr,max ~ m~^'^, 
significantly stronger than the {Nf = 2) universality prediction. In addition to that, if one 
calculates the chiral cumulants A(m, T) = Xct/x-k for the two different masses, they fail to 
cross in the regime of temperatures studied. 

This may be an indication that the quark mass m = 0.07A is already too small and one 
is entering the regime of "mesoscopic" QCD, where the mass is smaller or comparable to 
the smallest eigenvalue A^m = 0{1/V) of the Dirac operator. There is now a well developed 
theory of this regime p8| , |68[] , which we do not want to enter here. Instead, we have extended 



our calculations to larger masses ruu = rrid = (0.10, 0.13, 0.15, 0.17)A. In fig.l6 we show the 
maximum in the disconnected part of the scalar susceptibility as a function of the quark 
mass. For m > O.lOA the behavior is well described by a single exponent. Fitting the 
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dependence of the peak height on the hght quark mass, we find Xdis,max ~ ^ , quite 
consistent with the {Nf — 2) universahty prediction Xa,max ~ m~°'^^. 

VIII. CONCLUSIONS 

We have studied the structure of the instanton hquid at zero and finite temperature. 
For this purpose, we have introduced a method that allows us to extract the free energy of 
the instanton liquid from numerical simulations. Using this method, we can determine the 
instanton density self-consistently, by minimizing the free energy. 

We found that for the classical streamline interaction supplemented by a short range 
core, the instanton liquid in full QCD (with two light and one intermediate mass flavor) 
stabihzes at a density N/V — 0.174A^. Fixing the scale parameter in order to reproduce the 
"canonical" value N/V — Ifm^ (corresponding to (^G^) = (360 MeV)"^), we find a vacuum 
energy density e = — 280MeV/fm^ and a quark condensate (qq) = — (216MeV)^. The 
average instanton size is p = 0.42 fm, leading to a packing fraction / = 0.57r^p^/ {N/V) — 
0.14. 

This means that the interacting instanton ensemble, especially using the ratio ansatz 
required at finite temperature, is not as dilute as the random ensemble (where we fixed 
N/V = 1 fm~^ and p — 0.33 fm), which so successfully reproduces many hadronic correlation 
functions. The situation is quite analogous to the nuclear matter problem, where reproducing 
the saturation density and the binding energy requires a delicate balance between the short 
range repulsive core and the intermediate range attraction. In the instanton liquid, the 
strength of the core (the short range repulsion) controls the diluteness of the ensemble. Too 
much repulsion, however, leads to a very small total density of instantons, which can only 
be balanced by some intermediate range attraction. In the ratio ansatz, the short range 
core is weak and the intermediate range attraction fairly shallow, leading to a rather dense 
ensemble. 

We note that at this point, the instanton-antiinstanton interaction at short distances is 
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not a very well defined concept, and the core that is required to stabilize the ensemble (for the 
streamline ansatz) is a purely phenomenological parameter. In practice we have fixed this 
parameter in order to reproduce the instanton density and size distribution inferred from 



phenomenological considerations and observed in lattice simulations [|T7|,|T8|. Providing a 
more solid theoretical foundation for the instanton interaction will be an important direction 
for future study. One possibility, to derive the short range part of the instanton interaction 
from the cross section for isotropic multi-gluon production |Q, was already mentioned in 



section 4. Another idea, the inclusion of non-perturbative corrections to the beta function 



was discussed in |[T9[| . 



We have also extended our simulations of the instanton liquid to finite temperature. 
For this purpose, we have to replace the instanton interaction and fermionic overlap matrix 
elements with their finite temperature counterparts, satisfying the appropriate boundary 
conditions in a euclidean box with temporal extent 1/T. We have shown that in full QCD, 
there is a phase transition at T ~ 140 MeV in which the chiral symmetry is restored. The 
mechanism that drives the transition is the formation of polarized instanton-antiinstanton 
molecules. As a result, the number of small eigenvalues associated with unpaired topological 
charges is suppressed, and chiral symmetry is restored. Furthermore, for 2 light flavors the 
transition appears to be second order, while for Nf > 3 there is evidence for a weak first order 
transition. We have also investigated the phase diagram of QCD with many fermion flavors. 
Qualitatively, it is clear that many fermion flavors help the formation of molecules. As a 
result, the transition temperature drops as the number of light flavors is increased. For Nf 
larger than some critical value NJ"^, chiral symmetry is restored even at zero temperature. 
For Nc = 2,3 our results are consistent with A^|^** = Nc + C, where 1 < C < 2. We have 
also studied the location of the discontinuity of the quark condensate in the T — m plane. 

We have also determined the free energy associated with instantons, ignoring the pertur- 
bative (non-zero mode) contributions from quarks and gluons. At small T it rises slightly, 
but then starts to drop and becomes very small for T > 180 MeV. While the energy density 
and pressure above Tc are dominated by the perturbative Stefan-Boltzmann contributions, 
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instantons can provide a significant contribution to the so called "interaction measure" 
B = —{e — 3p)/4 even at T > Tg. Apart from the location of the phase transition, the most 
important result of our work is that B retains about half of its T = value at T = Tc. Using 
the trace anomaly, the interaction measure is related to the gluon and quark condensates. 
Indeed, analyzing lattice simulations of full QCD one can show that about half of the glue 



remains condensed at T = Tc ||69H71|. More specifically, lattice simulations show that the 
pressure remains small until T ^ 2Tc, while the energy density reaches the Stefan-Boltzmann 
value very quickly [^ , [72| , |57| . It was emphasized in that this behavior requires that only 



part of the bag pressure (or the gluon condensate) is removed across the phase transition. 
This is precisely the behavior observed in the instanton liquid model. 

Finally, we have studied the behavior of the Dirac spectrum and the mesonic suscepti- 
bilities near the phase transition. This is related to the question whether the behavior of 
these quantities is governed by critical scaling due to a nearby second order phase transition, 
and, more generally, what symmetries are restored in the transition. We find a peak in the 
scalar-isoscalar (cr) susceptibility, indicating the chiral phase transition. The dependence of 
the peak height on the light quark mass is consistent with the expected scaling behavior 
at fairly large quark masses m > O.IA, but is more singular at smaller masses, probably 
due to the onset of finite size effects. For m > O.IA, there is no hint of U{1)a restoration: 
there is no peak in the isovector scalar (5) susceptibility and the Dirac spectrum contains 
(approximate) zero modes above T^. At smaller masses the scalar-isovector (6) susceptibility 
shows an ehancement around the critical point, but with the available data it is difficult to 
say whether this behavior will persist for larger volumes. 
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APPENDIX A: INSTANTON INTERACTION AT ZERO TEMPERATURE 



In this appendix we specify the classical instanton interaction used in our simulations. 
From the form of the action for the two- instanton ansatz, one can easily show that the 
general form of the instanton-antiinstanton interaction in the case of the gauge group SU (2) 
is given by 

Sint = /9i(p) {so{R) + s,{R){u ■ Rf + S2{R){u ■ Rf) , (Al) 

where is the orientation vector introduced in section 2 and i?^ is the unit vector connecting 
the centers of the instanton and antiinstanton. The interaction is given in units of the 
single instanton action Sq = Pi{p)- The argument of the beta function is not uniquely 
determined without a higher order calculation. In practice we have taken the geometric 
average p = ^J pipA of the instanton radii. For large separation i?, the interaction ( |Al| ) can 
be embedded into SU (3) by making the replacement [u- Kf' \u- R\'^ and multiplying the 
first term by |m|2. 

For the ratio ansatz, the instanton-antiinstanton interaction can be parameterized by 



SiA ( 4.0 1.66 0.721og(r 



2^ 



(3i \ (r2 + 2.0)2 (1 + 1.68r2)3 (1 + 0.42r2) 

/ 16.0 2.73 \ 

^ I ~ (r2 + 2.0)2 ^ (1 + 0.33r2)3 ) ' 



+ TTT^^ihh l«r (A2) 



where r = R/ ^JpIPA is the instanton-antiinstanton separation in units of the geometric mean 
of their radii. In the ratio ansatz, the instanton-instanton interaction is non-zero if the two 
instantons have different color orientations. The interaction can be parameterized by 

Su 1 



/3i (1 + 0.43r2)3 

log(r2 



0.63|mP + 0.071|m|^ 



(A3) 

0.05|mP + 0.47|m| " 



(1 + 1.17r2)4 

In the streamline ansatz, conformal symmetry dictates that the interaction depends on the 
relative separation and the instanton radii only through the conformal parameter 
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The instanton-antiinstanton interaction is then given by 
SiA 4 



1/2 



(A4) 



4 1 - A* + 4A'log(A) 



Mp-4|M-i?|2 



(A5) 



A (A2 - 1)3 

+ 2 (l - + (1 + A^) log(A)) - 4|n ■ R\y + \u\^ + 2{uyiu*y 

As discussed in section 2, the instanton-instanton interaction vanishes in the streamhne 
ansatz. 



APPENDIX B: FERMIONIC OVERLAP MATRIX ELEMENTS AT ZERO 

TEMPERATURE 

A parameterization of the fermionic overlap matrix element between instantons and 
antiinstantons in the sum ansatz was already given in section 2 

TiA = tiu-R)^——^ -. (Bl) 

pipA (2.0 + R^/piPaY 

The matrix element in the ratio ansatz is very similar and we employ the same parameteri- 
zation. In the streamline ansatz, the matrix element depends on the separation and the radii 
again only through the conformal parameter A. The matrix element can be parameterized 
by IZl 

1 ciA^/^ 
TiA = i{u-R) — — — -— .^,3,^ , (B2) 

pipa (1 + i.25(A^ - 1) + C2(A^ - lyyi^ 

with ci = f and ca = 

APPENDIX C: INSTANTON INTERACTION AT FINITE TEMPERATURE 

In this appendix we give a parameterization of the classical instanton interaction at 
finite temperature. For reasons explained in section 4, there is no analog of the streamline 
interaction at finite temperature and we only specify the interaction in the ratio ansatz. 
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Up to small changes that have been introduced in order to improve the parameterization at 



small temperatures, the interaction is identical to the one given in |]52|. Here we also specify 
how to embed the SU{2) parameterization into SU{3). 

At finite temperature, the interaction is still at most quartic in the relative orientation 
vector u^. However, since four dimensional rotational invariance is broken, the interaction 
depends on |-U4p in addition to the invariants \u ■ -Rp and |mP appearing in the zero temper- 
ature interaction. For the same reason, the interaction can now depend separately on the 
spatial and temporal components of the vector i?^ = {R, R4) . At temperatures of interest 
here, the anisotropy in the dependence on i?^ turns out to be small. The parameterization 

SjA ^ 4.0 (3' _ ( 1.66 0.721og(r^) \ 

Pi (r2 + 2.0)2/32 + 5.21'''' 1^(1 + 1.68r2)3 + (l + 0.42r2)4j /32 + 0.75 ^ ' 

( 16.0 2.73 \ 
' + -, : — ^ttt —\u-RY' 



\ (r2 + 2.0)2 (l + 0.33r2)3y /?2^ 0.24 + 11.50r2/(l + 1.14r2) 

therefore depends only on r = R/^fpipA with R={R^ + Riy/\ The inverse temperature 
(3 = 1 / {T^piPa) is also given in units of the mean of the instanton radii. One can easily 
check that for T ^ the interaction ( PT|) reduces to the zero temperature ratio ansatz ( [A2D . 
The interaction between two instantons can be parameterized by 
Sii 1 (3^ 



/3i (1 + 0.43r2)3/52 + 5 33 

log(r2) (3^ 



0.63|Mr + 0.071|m 

2 I n /i'7l„-?|4 



(C2) 



;i + 1.17r2)4/52 + l.l7 



0.05|m|^ + 0.47|m| 



APPENDIX D: FERMIONIC OVERLAP MATRIX ELEMENTS AT FINITE 

TEMPERATURE 

At finite temperature, the fermionic overlap matrix element is still linear in the relative 
orientation vector m^. Due to the loss of Lorentz invariance at finite temperature, the matrix 
element depends separately on U4 and u ■ R 
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TiA = tUifi + I — ^ — /2. (Dl) 

A parameterization of the functions /i,2 was given in |5^. We have changed this parame- 
terization shghtly in order to improve the behavior at small temperatures. The result for /i 



IS 



(cosh(^)-cos(lf)+.;) 

1 

" ^) (e.p + -) + (I) (l - 0.69 exp {-IfL 



0.76/52 \ / /vrrX^ 0.18 



V (0.82r2 + 1)2(1 + 0.19/32)2^ \^ \PJ 1 + 0.12r2 

where r = \R\/ y/pJpA, r = R^/ ^fpTpA and /? = XjiT ^ pipA) are given in units of the mean 
instanton radius. The result for /2 is 

' vrr 







j cos 


(t) 


sinh 


cosh 1 




— cos 


f TXT 



+ K 



f^)(-p(-^) + f) + (l)(l + °-«-p(-^) 



4 = — Ktsft- m 



with Ki_2 given by 



One can easily verify that the parameterization ( |D 1[ - P4| ) reduces to ( [Bl| ) in the zero tem- 
perature limit. 



0.53 + ( CV 0.69 +(^ 



41 



REFERENCES 

[1] A. A. Belavin, A. M. Polyakov, A. S. Schwartz, Y. S. lyupkin, Phys. Lett. B59, 85 

(1975) 

[2] G. 't Hooft, Phys. Rev. 14D, 3432 (1976) 

[3] C. G. Callan, R. Dashen, D. J. Gross, Phys. Rev. D17, 2717 (1978) 
[4] E. V. Shuryak, Phys. Lett. B79, 135 (1978) 

[5] M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, Nucl. Phys. B 147 385, (1979) 
[6] E. V. Shuryak, Nucl. Phys. B203, 93, 116 (1982) 
[7] D. I. Diakonov, V. Yu. Petrov, NucL Phys. B245, 259 (1984) 
[8] E. M. Ilgenfritz, M. Miiller-Preussker, Nucl. Phys. B184, 443 (1981) 
[9] E. V. Shuryak, Nucl. Phys. B302, 559, 574, 599 (1988) 
[10] E. V. Shuryak, Nucl. Phys. B319, 521, 541 (1989) 

[11] E. V. Shuryak, J. J. M. Verbaarschot, Nucl. Phys. B410, 55 (1993); T. Schafer, 
E. V. Shuryak, J. J. M. Verbaarschot, Nucl. Phys. B412, 143 (1994) 

[12] T. Schafer, E. V. Shuryak, Phys. Rev. D50, 478 (1994) 

[13] T. Schafer, E. V. Shuryak, Phys. Rev. Lett. 75, 1707 (1995) 

[14] E. V. Shuryak, Rev. Mod. Phys. 65, 1 (1993) 

[15] M. C. Chu, J. M. Grandy, S. Huang, J. W. Negele, Phys. Rev. Lett. 70, 225 (1993); 
Phys. Rev. D 48, 3340 (1993) 

[16] B. Berg, Phys. Lett. B114, 475 (1981); M. Teper, Nucl. Phys. B20 (Proc. Suppl.), 159 
(1991) 

[17] M. C. Chu, J. M. Grandy S. Huang , J. W. Negele, Phys. Rev. D49, 6039 (1994) 

42 



[18] C. Michael, P. S. Spencer, Nucl. Phys. (Proc. Suppl.) B42, 261 (1995); Phys. Rev. D50, 
7570 (1995) 

[19] E. V. Shuryak, preprint. Stony Brook, SUNY-NTG-95-12 

[20] T. Schafer, E. V. Shuryak, preprint. Stony Brook, SUNY-NTG-95-23 

[21] E. V. Shuryak, J. J. M. Verbaarschot, Phys. Rev. D52, 295 (1995) 

[22] I. Zahed, Nucl Phys. B427, 561 (1994) 

[23] D. I. Diakonov, A. D. Mirhn, Phys. Lett. B203, 299 (1988) 

[24] M. A. Nowak, J. J. M. Verbaarschot, I. Zahed, Nucl. Phys. B325 (1989) 581. 

[25] R. D. Pisarski, L. G. Yaffe, Phys. Lett. B97, 110 (1980); D. J. Gross, R. D. Pisarski, 
L. G. Yaffe, Rev. Mod. Phys. 53, 43 (1981) 

[26] E.-M. Ilgenfritz, E. V. Shuryak, Nucl. Phys. B319, 511 (1989) 

[27] V. Khoze, A. Yung, Z. Phys. C50, 155 (1991) 

[28] E. V. Shuryak, M. Velkovsky Phys. Rev. D50, 3323 (1994) 

[29] M. C. Chu, S. Schramm, Phys. Rev. D51, 4580 (1995) 

[30] E.-M. Ilgenfritz, M. Miiller-PreuBker, E. Meggiolaro, Nucl. Phys. (Proc. Suppl.) B 42, 
496 (1995) 

[31] E.-M. Ilgenfritz, E. V. Shuryak, Phys. Lett. B325, 263 (1994) 

[32] T. Schafer, E. V. Shuryak, J. J. M. Verbaarschot Phys. Rev. D51, 1267 (1995) 

[33] T. Schafer, E. V. Shuryak, Phys. Lett. B356, 147 (1995) 

[34] A. V. Yung, Nucl. Phys. B297, 47 (1988) 

[35] J. J. M. Verbaarschot, Nucl. Phys. B362, 33 (1991) 



43 



[36] D. I. Diakonov, V. Yu. Petrov, Nucl. Phys. B272, 457 (1986) 

[37] M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, Nucl. Phys. B163 43, (1980) 

[38] H. Leutwyler, A. Smilga, Phys. Rev. D46, 5607 (1992) 

[39] E. V. Shuryak, J. J. M. Verbaarschot, Nucl. Phys. B341, 1 (1990). 

[40] E. V. Shuryak, O. V. Zhirov, Nuck. Phys. B242, 393 (1984) 

[41] J. Engels, J. Fingberg, F. Karsch, D. Miller, M. Weber, Phys. Lett. B252, 625 (1990); 
T. Blum, L. Karkkainnen, D. Toussaint, S. Gottlieb, Phys. Rev. D51, 5133 (1995) 

[42] E. V. Shuryak, Phys. Lett. B153, 162 (1985) 

[43] D. I. Diakonov, V. Yu. Petrov, JETP 89, 361 (1985); Phys. Lett. B147, 357 (1984) 
[44] D. I. Diakonov, V. Yu. Petrov, Nucl. Phys. B389, 109 (1991) 
[45] D. I. Diakonov, V. Yu. Petrov, Phys. Rev. D50, 266 (1994) 

[46] V. A. Novikov, M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, Nucl. Phys. B165, 67 

(1979) 

[47] M. A. Nowak, J. J. M. Verbaarschot, I. Zahed, Phys. Lett. B228, 115 (1989) 

[48] A. Smilga, J. Stern, Phys. Lett. B318, 531 (1993) 

[49] N. Dowrick, M. Teper, Nucl. Phys. (Proc. Suppl.) B42, 237 (1995) 

[50] J. J. M. Verbaarschot, Act. Phys. Pol. 25, 133 (1994); Nucl. Phys. B427, 534 (1994) 

[51] B. J. Harrington, H. K. Shepard, Phys. Rev. D17, 2122 (1978); B. Grossman, 
Phys. Lett. A61, 86 (1977) 

[52] E. V. Shuryak, J. J. M. Verbaarschot, Nucl. Phys. B364, 255 (1991) 

[53] R. D. Pisarski, F. Wilczek, Phys. Rev. D29, 338 (1984) 



44 



[54] F. Wilczek, Int. J. Mod. Phys. A7, 3911 (1992); K. Rajagopal, F. Wilczek, 
Nucl. Phys. B399, 395 (1993) 

[55] A. Kocic, J. Kogut, Phys. Rev. Lett. 74, 3109 (1995) 

[56] F. Karsch, E. Laermann, Phys. Rev. D50, 6954 (1994) 

[57] C. DeTar, Quark Gluon Plasma in numerical simulations of lattice QCD, to appear in 
"Quark Gluon Plasma 2", R. Hwa, ed., World Scientific (1995) 

[58] T. Banks, A. Zaks, Nucl. Phys. B82, 196 (1982) 

[59] E. V. Shuryak, Phys. Lett. B193, 319 (1987) 

[60] Y. Iwasaki, K. Kanaya, S. Sakai, T. Yoshie, Nucl. Phys. (Proc. Suppl.) B34, 314 (1994) 
[61] N. Christ, Nucl. Phys. (Proc. Suppl.) B30, 323 (1993) 
[62] N. Seiberg, Phys. Rev. D49, 6857 (1994) 

[63] S. Chandrasekharan, Nucl. Phys. (Proc. Suppl.) B42, 475 (1995) 

[64] J. Engels, F. Karsch, H. Satz, Nucl. Phys. B315, 419 (1989) 

[65] H. Leutwyler, Phys. Lett. B284, 106 (1992) 

[66] T. Banks, A. Casher, Nucl. Phys. B169, 103 (1980) 

[67] E. V. Shuryak, Comm. Nucl. Part. Phys. 21, 235 (1994) 

[68] J. J. M. Verbaarschot, 1. Zahed, Phys. Rev. Lett. 70, 3852 (1993); J. J. M. Verbaarschot, 
Phys. Rev. Lett. 72, 2531 (1994) 

[69] V. Koch, G. E. Brown, Nucl. Phys. A560, 345 (1993) 

[70] Y. Deng, Nucl. Phys. (Proc. Suppl.) B9, 334 (1989) 

[71] C. Adami, T. Hatsuda, I. Zahed, Phys. Rev. D43, 921 (1991) 



45 



[72] J. B. Kogut, D. K. Sinclair, K. C. Wang, Phys. Lett. B263, 101 (1991) 
[73] E. V. Shuryak, J. J. M. Verbaarschot, Phys. Rev. Lett. 68, 2576 (1992) 



46 



FIGURES 

FIG. 1. Classical instanton-antiinstanton interaction in the streamline (dash-dotted line) and 
ratio ansatz (short dahed line). The interaction is given in units of the single instantanton action 
So for the most attractive (cos^ = 1) and most repulsive (cos^ = 0) orientations. The dash-dotted 
curves show the original streamline interaction, while the solid curves show the interaction includ- 
ing the core introduced in section 3. Fig. lb shows the fcrmionic overlap matrix elements in the 
streamline (solid curve) and ratio ansatz (dashed curve). The matrix elements are given in units 
of geometric mean of the instanton radii. 

FIG. 2. Free energy, average instanton size and quark condensate as a function of the instanton 
density in the pure gauge theory. All quantities are given in units of the scale parameter Aqcd- 

FIG. 3. Free energy, average instanton size and quark condensate as a function of the instanton 
density in the theory with two light and one intermediate mass flavor. All quantities are given in 
units of the scale parameter Aqcd- 

FIG. 4. Distribution of instanton sizes, eigenvalues spectrum of the Dirac operator and distri- 
bution of fermionic overlap matrix elements in pure gauge theory, p, X and Tja are given in units 
of the scale parameter. The distribution functions have arbitrary units. 

FIG. 5. Distribution of instanton sizes, eigenvalue spectrum of the Dirac operator and distri- 
bution of fermionic overlap matrix elements in full QCD. 

FIG. 6. Trajectory of the quark condensate (in units of fm~^ for Nf=2 in the critical region, 
T = 150 MeV. We also show the corresponding spectrum of Dirac eigenvalues A [fm~^]. The 
normalization of the eigenvalue distribution is arbitrary. 

FIG. 7. Double logarithmic plot of the quark condensate (in units fm~^) versus the valence 
quark mass for different number of light flavors (with dynamical mass m = 20 MeV). Fig.7b 
shows the same plot for Nf = 5 and different quark masses m. 
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FIG. 8. Schematic phase diagram of the instanton Hquid for different numbers of quark flavors, 
Nf=2,3 and 5. We show the phase of chiral symmetry breaking in the temperature-quark mass 
planes. For Nf=2, open squares indicate points where we found large fluctuations of the chiral 
condensate, the cross indicates the approximate location of the singularity. In the two other figures 
the open squares correspond to points where we find a plateau in the valence mass dependence of 
the chiral condensate, while solid squares correspond to points where such a plateau is absent. The 
crosses and the dashed lines connecting them show the approximate location of the discontinuity 
line. 

FIG. 9. Instanton density, free energy and quark condensate as a function of the temperature 
in full QCD with two light and one intermediate mass flavor. 

FIG. 10. Typical instanton ensembles for T = 75, 123 and 158 MeV. The plots show projec- 
tions of a four dimensional (3A~^)^ x box into the 3-4 [z axis-imaginary temperature) plane. 
Instantons and antiinstanton positions are indicated by -|- and — symbols. Dashed, solid and thick 
solid lines correspond to fermionic overlap matrix elements Tja > 0.40, 0.56, 0.64, respectively. 

FIG. 11. Trajectories of the quark condensate at three different temperatures T = 75, 130, 158 
MeV. The quark condensate is given in units of the scale parameter. 

FIG. 12. Spectrum of the Dirac operator for different temperatures T = 75, 130, 158 MeV. 
The normalization of the spectra is arbitrary (but identical) . 

FIG. 13. Distribution of the color orientation angle cos^ a = |«4p/|'Up for different tem- 
peratures T = 75, 130, 158 MeV. The distributions are normalized to the corresponding SU{3) 
measure. The overall scale is arbitrary. 
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FIG. 14. Mesonic susceptibilities as a function of temperature for quark masses 
nT-u = nT-d = O.lOA and rus = 0.70A. The scalar a, scalar 6 and pseudoscalar tt susceptibilities 
are denoted by solid circles, solid squares and open squares, respectively. The lower panel shows 
the disconnected part of the scalar susceptibility and the difference Xn — XS- The solid lines show 
a smooth fit to the data. 

FIG. 15. Mesonic susceptibilities for a smaller light quark mass = rud = 0.07A. Curves 
labeled as in figure 14. 

FIG. 16. Maximum of the disconnected part of the scalar susceptibility as a function of the 
quark mass. Also shown is the parameterization x ~ m~°'^^. 
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TABLES 

TABLE I. Parameters of the pure gauge instanton ensemble in the streamline ansatz for various 

vahics of the core parameter A. The density and average instanton size are given in units of Aqqd 
and fm. The dimensionless scale parameter / is introduced in the text. 
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